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The structure, magnetic properties, and lattice dynamics of ordered Fe-Pt alloys with three stoi¬ 
chiometric compositions, FeaPt, FePt and FePt 3 , have been investigated using the density functional 
theory. Additionally, the existing experimental data have been complemented by new measurements 
of the Fe projected phonon density of states performed for the FeaPt and FePt 3 thin films using 
the nuclear inelastic scattering technique. The calculated phonon dispersion relations and phonon 
density of states have been compared with the experimental data. The dispersion curves are very 
well reproduced by the calculations, although, the softening of the transversal acoustic mode TAi 
leads to some discrepancy between the theory and experiment in FeaPt. A very goood agreement 
between the measured spectra and calculations performed for the tetragonal structure derived from 
the soft mode may signal that the tetragonal phase with the space group PA/mbm plays an im¬ 
portant role in the martensitic transformation observed in FeaPt. For FePt 3 , the antiferromagnetic 
order appearing with decreasing temperature has been also investigated. The studies showed that 
the phonon density of states of FePt 3 very weakly depends on the magnetic configuration. 


I. INTRODUCTION 

The Fe-Pt alloys constitute a very important class 
of materials because of their various interesting physi¬ 
cal properties and high application potential. Crystal- 
lographically, they can be either chemically disordered 
fee or ordered fcc/fct structures. According to the Fe-Pt 
phase diagranfS, at high temperatures an fee solid so¬ 
lution of the components is observed. With decreasing 
temperature this disordered Al-type structure exhibits 
order-disorder transformations leading to three different 
ordered phases with structures and physical properties 
depending on the chemical composition. 

At temperatures below 1570 K, the tetragonal Llo 
structure is formed for the almost equiatomic concentra¬ 
tion region from approximately 35 to 55 atomic percent of 
Pt. In Fe-Pt alloys with lower or larger Pt concentration, 
the formation of the stable cubic LI 2 structures, FePt 3 
and Fe 3 Pt, is observed at temperatures below 1120 K 
and 1620 K, respectively. Moreover, in the region below 
670 K the Fe-rich alloys show the anomalously low ther¬ 
mal expansion coefficient (Invar effect)^ and undergo a 
martensitic transformation . 3 The martensite struture de¬ 
pends on the Fe concentration and the chemical order. In 
the ordered state, the Ll 2 Fe 3 Pt phase is stable down to 
approximately 60 K, whereas the disordered FeaPt starts 
to transform to a bcc martensite already at room tem¬ 
perature. As a precursor of the distorted phase, the soft¬ 
ening of the transverse acoustic phonons is observed in 
the [ 110 ] direction , 4 however, the role of the soft mode 
in the mechanism of the structural transition is not fully 
explained. 

The ordered equiatomic FePt alloy with the Llo struc¬ 
ture is composed of alternately stacked layers of Fe and 
Pt atoms along c-axis (AuCu-type structure with the 
PA/mmm space group). This structure exhibits ferro¬ 
magnetic (FM) order below T c = 750 K. The FePt alloys 
with stoichiometry around 1:3 can form a cubic phase 


FeaPt or FePt 3 (with the Pm3m symmetry). In FeaPt 
(or FePta), the Fe (Pt) atoms occupy the cube corners 
and the Pt (Fe) atoms occupy the face-center positions. 

In Fe-Pt alloys, both Fe and Pt atoms carry mag¬ 
netic moments, however, the Fe moments are significantly 
larger. The orientation of local magnetic moments de¬ 
pends on the Fe and Pt concentrations, arrangement of 
alloy components and temperature. The disordered Al 
phase of FePt and FePta is ferromagnetically ordered in 
contrast to the Fe 3 Pt disordered alloy that is paramag¬ 
netic. Differently, the ordered FePta alloy is paramag¬ 
netic while FePt and Fe 3 Pt are ferromagnetic. Moreover, 
below 170 K, FePt 3 is antiferromagnetic (AFM-I phase) 
with wavevector qi = — (|,|,0)P Subsequently, below 
80 K it undergoes a second phase transition into an AFM- 
II phase with wavevector (72 = — (^,0, 0)Plt has been 
discussed,^ that the latter AFM state is metastable and is 
probably induced by antiphase boundaries or small com¬ 
positional variations. Such a wide variation of magnetic 
structures in the FePt alloys is evidently a consequence of 
different atomic configurations around Fe atoms, which 
in turn, have a considerable effect on the electronic struc¬ 
ture of these alloys. 

Differences in crystal structures and chemical composi¬ 
tions influence also lattice dynamics of Fe-Pt alloys. Ex¬ 
perimentally, the knowledge on phonon spectra for Fe-Pt 
alloys comes from fitting the phonon dispersion relations 
obtained from the inelastic neutron scattering (INSjMal 
to the Born-von Karman model and from nuclear inelas¬ 
tic scattering (NIS) measurement:!®®. The latter tech¬ 
nique provides information on the Fe contribution to the 
phonon density of states (PDOS), and it has been ap plied 
to ordered and disorder polycrystalline Fe-P t alloys}®® 
thin films!® as well as to nanocluster!! 3 ® The theo¬ 
retical studies based on ab initio methods enabled cal¬ 
culation of the phonon dispersion relations of Fe/Pt 
multilayer!®, estimation of the magnetic anisotropy of 
Fe/Pt(001) and Pt/Fe/Pt(001) system!®, identification 
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of the structural and magnetic phases in FePt surface 
alloyfi®- ; , determination of the partial Fe PDOS in or¬ 
dered FePlP-^ and investigation of the soft mode behav¬ 
ior in Fe 3 PtP-^ For FePt, the first-principles results were 
compar ed wit h the partial Fe PDOS measured by the NIS 
method. 14 ^Nevertheless, no direct comparison between 
the calculated and measured dispersion curves has been 
performed for any of the Fe-Pt alloys. 

In this paper, we present the results of ab initio studies 
on the structural, magnetic, and dynamical properties for 
three ordered alloys: FePt 3 , FePt, and FeaPt. The ob¬ 
tained results are confronted with previously measured 
phonon dispersion curves and PDOSs as well as new NIS 
measurements performed on FeaPt and FePta thin films. 
Using the same formalism and software for all these sys¬ 
tems allows us to make detailed comparison and discuss 
the differences in lattice dynamics of different Fe-Pt al¬ 
loys that otherwise would not be possible. In particular, 
we are studying the influence of local atomic configura¬ 
tions and magnetic order on phonon spectra. For FeaPt, 
we analyse the low-symmetry phase obtained from the 
soft mode, and compare the calculated phonon spectra 
with the low-temperature measurements. 

The paper is organized as follows. In Section [IlJ the 
structure, magnetic ordering, and dynamical properties 
of room temperature phases of ordered Fe-Pt alloys are 
analysed and compared with the experimental data. In 
Section |III[ the detailed analysis of the anisotropy of 
phonon density of states and the low-temperature phases 
of FePta and FeaPt is presented. Section IV concludes 
the results. 


II. CALCULATIONS VS MEASUREMENTS 
A. Structure and magnetic ordering 

Three ordered phases of Fe-Pt alloys were modelled 
by imposing the symmetry restrictions of the P4/mmm 
(Ll 0 ) and Pm3m (Ll 2 ) space groups on the crystal struc¬ 
ture. Some results for the Llo equiatomic structure have 
already been presented in our previous paperf 1 ^ The cal¬ 
culations of cubic phases have been performed using the 
same technique. Structure optimization was achieved 
using the VASP package.^ The spin-polarized density 
functional theory (DFT) calculations were carried out 
within the generalized gradient approximation us ing the 
Perdew, Burke, and Ernzerhof (PBE) functional “ The 
wave functions were sampled according to Monkhorst- 
Pack scheme with a k -point mesh of (4,4,4). The struc¬ 
tural calculations were performed on a 2 x 2 x 2 super¬ 
cell (containing 32 atoms) with periodic boundary con¬ 
ditions. During the structure optimization only lattice 
constants are modified, none of the atoms change their 
position in the unit cell as all of them are placed in high 
symmetry crystallographic sites. 

The calculated lattice parameters are presented in Ta¬ 
ble I together with the Strukturbericht and space group 


TABLE I. The Strukturbericht symbols, crystallographic 
space groups, and lattice parameters of the ordered Fe-Pt 
alloys with stoichiometric concentrations. The experimental 
values of lattice constants are presented in parenthesis. 



Struktur¬ 

bericht 

symbol 

Space 

Group 

symbol 

Lattice 

parameters 

(A) 

FePt (FM) 

Llo 

P4/mmm 

a = 3.838 (3.852) 
c = 3.739 (3.713) 

Fe 3 Pt (FM) 

Ll2 

Pm3m 

a = 3.737 (3.750) 

FePt 3 (FM) 

Ll2 

Pm3m 

a = 3.914 (3.866) 

FePt 3 (AFM) 

Ll2 

Pm3m 

a = 3.891 (3.866) 

FePt 3 (NM) 

Ll2 

Pm3m 

a = 3.877 (3.866) 


TABLE II. The magnitude of the magnetic moments in [ib 
calculated for three ordered Fe-Pt alloys with stoichiometric 
concentrations. 


Atom 


Fe 3 Pt (FM) FePt (FM) FePt 3 (FM) FePt 3 (AFM) 


Fe 

Pt 


2.711 

0.349 


2.966 

0.354 


3.258 

0.359 


3.325 

0.000 


symbols describing these structures. The determined val¬ 
ues differ slightly from experimental data placed in the 
parenthesis. The FePt and FeaPt structures have been 
calculated assuming ferromagnetic arrangement of the lo¬ 
cal magnetic moments. This is not the case for FePta 
which is paramagnetic at room temperature. To investi¬ 
gate the influence of magnetic ordering on the lattice pa¬ 
rameter of FePta, we have performed calculations for the 
FM, AFM, and nonmagnetic (NM) phases. Simple AFM 
ordering results in lattice parameter reduction, in com¬ 
parison with the FM arrangement, however, the shortest 
lattice constant for FePta structure has been obtained 
within the NM calculations. 

The magnetic moments calculated for the ordered Fe¬ 
Pt structures observed at room temperature are pre¬ 
sented in Table II. With increased iron concentration, 
magnetic moments on both types of atoms decrease and 
the calculated values are in goo d agreement with the re¬ 
sults presented previously! 22 * 23 ' The structure and mag¬ 
netic ordering of the low-temperature phases are dis¬ 
cussed in further sections. 


B. Phonon dispersion relations 

The phonon dispersion relations and PDOS were cal¬ 
culated with the direct method implemented in the 
PHONON codeP^ This method utilizes the Hellmann- 
Feynman (H-F) forces obtained by performing small 
atomic displacements of nonequivalent atoms from their 
equilibrium positions. From them the dynamical matrix 
is determined and diagonalized to obtain the phonon fre¬ 
quencies at each wave vector. Since the DFT calculations 
are related to T = 0 K, we have adopted the experimen- 
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FIG. 1. (Color online) The phonon dispersions along the high 
symmetry directions calculated for the Llo FePt structure. 
The experimental data are taken from Ref. [6|. 


tal lattice constants at room temperature for comparison 
between the measured and calculated phonon frequen¬ 
cies. 

The experimental dispersion curves of the Fe-Pt al¬ 
loys have been measured using INS and published some 
decades agoP^ The measurements were done for the rel¬ 
evant monocrystals with a high chemical order at room 
temperature. The phonon frequencies have been mea¬ 
sured along the high symmetry directions of the cubic 
symmetry: [1 0 0], [1 1 0] and [1 1 1], and additionally in 
[0 0 1] and [1 0 1] in the Ll 0 tetragonal phase. In the Ll 2 
cubic phase, the latter two directions are equivalent to [1 
0 0] and [110], respectively. For FesPt, the temperature 
and concentration dep ende nce of the soft mode were also 
studied experimentally.^^ 

The conventional cell of the ordered Llo FePt structure 
contains four atoms, but it can be reduced to a prim¬ 
itive tetragonal cell with unchanged c lattice constant 
and a constant reduced to a^/2/2. The primitive cell 
contains two atoms only, thus 6 branches of phonon dis¬ 
persion curves are expected. They are presented in FigjT] 
in comparison with the data measured by INS along the 
main crystallographic directions of the conventional cell 
(taken from Ref. |6j . The agreement between measured 
and calculated values is very good. The most signifi¬ 
cant discrepancies are observed for the high-frequency 
branches. For example, the calculated high-frequency 
dispersion curves along (£,£,0) and (£,£,£) directions are 
less dispersive than the related experimental curves. 

The primitive cell of the Ll 2 structure of Fe 3 Pt or 
FePt 3 contains 4 atoms. Hence, one can expect 12 
phonon dispersion curves. They are shown in Fig. [2] Dif¬ 
ferent symbols of the experimental points correspond to 
two types of phonon polarizations: dark for longitudinal 
and light for transversal vibrations. In the top panels^ the 
experimental data taken for the FePt 3 alloy at 295 KO are 
compared with dispersion relations calculated for three 
different magnetic phases: FM, AFM, and NM states. 
The frequencies of the AFM phase are slightly higher 
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FIG. 2. (Color online) The phonon dispersions along the high 
symmetry directions calculated for the Ll 2 structures of FeaPt 
and FePt 3 . The experimental data are taken from Refs. [?] 
and 2]. 


than those of the FM phase but differences between the 
dispersion curves are negligible. Both of them reproduce 
the experimental data quite well, especially the frequen¬ 
cies calculated at the T point are in very good agreement 
with measured values. The most significant discrepan¬ 
cies are observed for the highest optical curves close to 
the X and R points. The NM calculations improve the 
agreement in the high-frequency region close to the X 
point, but generally the NM state generates larger dis¬ 
crepancies than the FM or AFM states. Each kind of 
calculations leads to too high frequencies at the R point. 
Probably, the experimental signals were not registered 
there because of too low intensity in the Brillouin zone 
chosen in the experiment. 

Moreover, in the measured curves some steep wiggles 
of the recorded frequencies of the longitudinal acoustic 
and the neighboring optical branches along the (0,0,£) 
direction have been noticed 7 ! For that acoustic mode, 
at the X point, our calculations generate a frequency of 
4 THz which is higher than 3.4 THz suggested from the 
fitting procedure of the Born-von Karman model. 7 Addi¬ 
tional analysis of the transverse acoustic branch along the 
(£,£,£) direction have been made as the downward bend 
of the curve measured close to the R point could indicate 
the possible tetragonal instability at low temperatures.^ 
Our calculations do not show any anomalous behavior at 
the R point that could confirm this suggestion. 
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In the calculations performed for the Fe-rich LI 2 struc¬ 
ture, the imaginary frequencies appear for the transver¬ 
sal acoustic mode along the (£,£, 0 ) direction and ap¬ 
proach the value 2.8i THz at the M point (bottom panel 
of Fig{2]). This behavior agrees with the measurements 
taken at 440, 295 and 80 K that showed significant fre¬ 
quency softening of the transversal acoustic (TA 4 ) mode 
registered along the (£,£, 0 ) direction with decreasing 
temperatureP In the bottom panel of Fig{2j the disper¬ 
sion curves measured at 295 K (dark and light points) are 
shown together with the TA-| branch measured at 295 K 
and 80 K (red circles and diamonds, respectively). The 
soft mode is described by the M 2 representation and is 
related to the vibration of Fe atoms onlyP In Ref. 8 j 
the authors noticed also the other mode (M 4 ) with the 
Fe atoms vibrating in a similar manner. The measured 
frequency of that optical mode is equal to the frequency 
of the M 2 mode at temperature of 440 K, and it stays 
almost constant with decreasing temperatureP In the 
present calculations, the M 4 mode vibrates with consid¬ 
erably higher frequency and this result confirms previous 

calculations Pi! 

In calculations, the low-frequency acoustic branch 
along the (£,£,£) direction slightly softens in the vicin¬ 
ity of the r point. The other acoustic branches and the 
high-frequency optical branches (between 6 and 8 THz) 
are well reproduced. However, the position of the low- 
frequency optic branches (close to 4 THz) is intriguing. 
In this frequency range, all calculated optic branches are 
located significantly lower than the experimental data. 
Even at the T point, almost 1 THz discrepancy is ob¬ 
served. Calculation errors can be excluded since the 
phonon dispersion relations are in very good agreement 
with the curves calculated previously by Gruner et all 19 1 
The shift of phonon branches seems to be related to the 
martensitic phase transformation as well as the imagi¬ 
nary frequency at the M point. The condensation of the 
soft mode drives the structure to a lower symmetry which 
is described in Sec. ITTTcl 


C. Phonon density of states 

In previous studies, the phonon density of states has 
been deduced from the measured dispersion relations us¬ 
ing the harmonic approximation and the Born-von Kar- 
man modelP? The fitting procedure of the model to the 
measured data enables to determine the atomic force con¬ 
stants. Using the appropriate summations the total or 
partial (for different kinds of atoms) PDOS have been 
obtained for the Ll 0 FePt structure.^ The vibrational 
frequencies cover the region from 2 to 8 THz (8 to 32 
meV) and as it can be expected the contribution of the 
heavier atoms to the total PDOS is dominant at low fre¬ 
quencies and that of the lighter atoms is important at 
high frequencies. 

The total PDOS of FesPt and FePt 3 have been shown 
and discussed in Ref. [81 The main difference between 
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FIG. 3. (Color online) The Fe partial PDOS of FeaPt and 
FePt 3 measured by NIS compared with the calculated spectra. 
The experimental data are taken from the NIS measurements 
performed as a part of this work. 


them is the high frequency limit of the spectra that is 6 
and 8 THz for FePt 3 and FesPt, respectively. The partial 
spectra of those structures have not been calculated from 
the measured dispersion curves. 

The information on vibrational frequencies in Fe¬ 
Pt alloys h as also been obtained from the NIS 
measurements P-Q 12 b 5 I For FesPt, the temperature depen¬ 
dent spectra were collected to study the low-temperature 
anomalies in that structure ? 1 For FePt 3 , the NIS spectra 
of 57 Fe atoms were measured on crystalline alloys that 
were chemically disordered, partially ordered, and Ll 2 
ordered . 12 First experiments done on the bulk crystalline 
FesPt and FePt 3 samples provided the PDOS averaged 
over all crystallographic directions ? 1 ^The development 
of NIS technique accompanied with the layer by layer epi¬ 
taxial preparation method al lows t o record the data even 
on monolayers of 57 Fe atoms? 2 ® 23 An appropriate choice 
of a substrate of the FePt thin films ordered in the Llo 
structure facilitated to perform the measurements of Fe 
partial PDOS along the (1,0,0) and (1,0,1) directions.^ 
These experiments have shown the asymmetry of spectra 
in good agreement with the first principles calculations. 

Currently, to complete the experimental data, the Fe 
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partial PDOS of FeaPt and FePt 3 have been measured 
by NIS at the P01 beamline of Petra III (Hamburg, Ger¬ 
many). The samples consisted of 30 nm films grown on 
MgO(lOO) substrates by co-evaporation of 57 Fe (98% iso¬ 
topic enrichment) and Pt at a temperature of 500°C. The 
measurements have been carried out at room tempera¬ 
ture in grazing incidence geometry. The resolution of the 
monochromator that defines the actual energy resolution 
of the measured data was 1 meV. 

In Fig. [3j the measured Fe projected PDOSs are shown 
together with the data obtained from the calculations. In 
the case of the FeaPt structure, the experimental phonon 
spectrum covers the wide frequency range from 2 up to 
9 THz. The disagreement between the experimental and 
calculated data is mainly observed in the low-frequency 
region (2-4 THz), where the contribution of calculated 
PDOS is much larger than experimental one. This ef¬ 
fect is related to the softening of modes calculated for 
the high-temperature austenite phase at 0 K. The high 
frequency end of the FeaPt vibrational spectrum is well 
reproduced. In the measured spectrum, a clearly vis¬ 
ible shoulder at ~2 THz corresponds to the TAa soft 
mode identified previously in the FeaPt bulk sampld 8 and 
nanoclusters 

In FePta, Fe atoms contribute to the PDOS in a much 
narrower range of frequencies with the main band centred 
around 6 THz, which is very well reproduced in the cal¬ 
culations. Four sharp peaks noticeable in the theoretical 
spectrum are weakly visible in the experiment because of 
the broadening of each phonon line. The experimental 
resolution cannot explain such large broadening that is 
rather caused by the anharmonicity of a material or size 
and disorder effects in a measured specimen. The cal¬ 
culations are based on the harmonic approximation and 
performed for a perfect bulk crystalline sample. This 
implies an infinite phonon lifetime, which is not a likely 
situation to be found in a nanometric size sample of Fe- 
Pt alloys. The sample size, grain boundaries, and crystal 
defects reduce the lifetime. The effect of finite phonon 
lifetime can be approximated by applying a damped har¬ 
monic oscillator (DHO) model to the ab initio calculated 
phonon spectrum 12 . This phenomenological model de¬ 
scribes the ’’damping” of the phonons i.e. it provides 
a finite phonon lifetime, which leads to broadening of 
phonon lines. Mathematically, this process is described 
as a convolution of the original PDOS with an appropri¬ 
ate broadening function D(E,E*): 


F DHO (E*) = J D{E,E*)F(E)dE (1) 


It has been shown that the damped harmonic oscillator 
function is a good choice for D(E, E*)&^. 


D(E, E*) 


1 1 


( 2 ) 


where Q is the quality factor, which describes the num¬ 
ber of periods a phonon will remain in a given state. 


Hence, the lower the value of Q, the shorter the lifetime 
and the larger broadening. This function is essentially a 
Lorentzian with a width that increases with E*, mean¬ 
ing that the PDOS features become broader at higher 
energy. The final PDOS (dashed blue line in Fig. [3]) is 
thus constructed by convolution of the theoretical PDOS 
with the DHO function using Eq. §■ In this model, only 
the Q factor is adjusted to provide the best match be¬ 
tween the theoretical and experimental PDOS. We fitted 
the calculation to the experimental data for FeaPt using 
Q=30, which provides a better match of the high energy 
peak intensities. This rather large value of Q indicates 
that only little damping takes place in that particular 
thin film system. However, this broadening does not im¬ 
prove the agreement between the experiment and theory 
at lower frequencies that is caused by the soft-mode be¬ 
havior. For FePta, a very good agreement is observed for 
Q=12, which indicates significant reduction of phonon 
lifetime in this system. 


III. DISCUSSION 
A. Anisotropy of PDOS 

Using first principles calculations one can determine 
the total PDOS as well as the partial PDOS for different 
atoms and directions. In Figs. [4]and[5j the PDOS calcu¬ 
lated for FePt, FeaPt and FePta are shown. Besides the 
total PDOS, the partial spectra of each type of atom pro¬ 
jected on the three orthogonal directions [100], [010] and 
[001] are calculated. Due to atomic site symmetries some 
spectra are equivalent. For the Llo tetragonal structure, 
the spectra taken along [100] and [010] are equivalent 
and their sum calculated separately for Fe or Pt atoms 
is shown in Fig. [4] as Fe X!/ or Pt xy , respectively. Sim¬ 
ilarly, the spectra projected on the [001] direction are 
named Fe 2 or Pt 2 . The strong asymmetry of Fe vibra¬ 
tions shown in Fig . [4] has already been observed in the 
NIS experiment P21 The vibrations along two directions 
are separated almost completely and the vibrations in the 
z direction contribute to a high frequency peak around 
7.2 THz. In contrast, the Pt atoms vibrations along the 
z axis dominate in a narrow frequency range from 2 to 3 
THz, while those in the xy-plane are spread in the whole 
spectrum including the peak at the lowest frequency. 

For the LI 2 structures, there are two kinds of atomic 
sites, atoms placed at corners with cubic site symme¬ 
try and atoms situated in the center of cubic faces with 
tetragonal symmetry. For Pt (or Fe) atoms placed in 
cubic symmetry sites of the FeaPt (or FePta) unit cell 
the PDOS projected on three crystallographic directions 
are equivalent. For the remaining three Pt (or Fe) atoms 
placed in the center of a cubic face, two equivalent direc¬ 
tions lying in each unit face are discriminated as ”in Pt 
(or Fe) plane” and the third direction, perpendicular to 
the unit face is called ’’out of Pt (or Fe) plane” (Fig§- 
The negative values of frequencies shown in the left pan- 
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FIG. 4. (Color online) The phonon density of states calculated 
for the Llo FePt structure. 

els of Fig[5]are related to the soft mode existing in Fe 3 Pt. 
The soft mode dominantly consists of Fe vibrations in the 
’’out of Pt plane” direction. The Pt atom vibrations are 
mainly observed in the narrow frequency range between 2 
and 4 THz. There is a significant difference between the 
Fe partial PDOS collected along the ”in Pt plane” and 
’’out of Pt plane” directions. The ’’out of plane” vibra¬ 
tions are involved in the low frequency region and the ” in 
plane” vibrations are distributed at higher frequencies. 

The spectrum calculated for Fe atoms in FePt 3 is char¬ 
acterized by a high frequency peak at about 6 THz with 
high intensity and two significantly smaller peaks at 4 
and 5 THz. The pronounced differences are observed for 
the PDOS calculated for Pt atoms vibrating along the 
direction parallel (”in Fe plane”) or perpendicular to the 
cell face (’’out of Fe plane”). 

Comparing the spectra for the three Fe concentrations 
one can see that the highest frequency of the total PDOS 
decreases with decreasing Fe content. 

B. Low-temperature AFM phases of FePt 3 

The chemically disordered phase of the Fe-Pt alloy 
with stoichiometry 1:3 is ferromagnetic with Fe moments 


of 2 /is. With decreasing temperature, the alloy trans¬ 
forms from the disordered fee to ordered LI 2 structure 
with paramagnetic properties. Further decreasing the 
temperature leads to the arrangement of magnetic mo¬ 
ments within two different types of AFM order. In the 
bulk, the appearance of these two AFM phases dep ends 
strongly on the deviation from the ideal composition! 27 ! 28 ! 
The FePt 3 phase is maintained in Fea,Pti_ x alloys for 
0.22<x<0.4lP2l Upon cooling, the AFM-I order devel¬ 
ops below 160 K, where the Fe moments order in al¬ 
ternating ferromagnetic (110) planes. This transition is 
of second order. The AFM-I tetragonal unit cell is ob¬ 
tained by doubling the cubic unit cell with the basic vec¬ 
tors at=(l,l,0)a c , 6t=(l,-l,0)a c and Ct=(0,0,l)a c . The 
moments carried by the Fe atoms are equal to 3.3 /x b■ 
With decreasing temperature, this phase transforms sub¬ 
sequently to the AFM-II phase where the ferromagnet- 
ically ordered (001) planes are stacked antiferromagnet- 
ically. It results in doubled lattice constant in the c- 
direction and the orientation of tetragonal unit cell is 
a t =(l,0,0)a c , 5 t =(0,l,0)a c and c t = (0,0,2)a c . The sym¬ 
metry of both phases is characterized by the P4/mmm 
space group, however, the magnetic unit cells and posi¬ 
tions of non-equivalent atoms are different (Table III). 

In Table III, the structural data, magnetic moments 
and calculated energy of the AFM-I and AFM-II super¬ 
cells are compared with the artificial FM, AFM, and NM 
structures of FePt 3 . The calculations performed for both 
AFM phases show that the ground state energy of the 
AFM-I phase is slightly lower than that of AFM-II phase, 
however, the energy difference between AFM-I and AFM- 
II is not appreciable. Similarly, the energies of artificial 
high-temperature structures of FM and AFM magnetic 
order are almost the same. The magnitudes of magnetic 
moments of Fe atoms vary from 3.258 /xb to 3.325 /is 
for the cubic phases and from 3.272 /xb to 3.306 /xb for 
tetragonal AFM structures. As it is shown in Table III, 
the lattice parameters of FePt 3 in AFM-I stay almost un¬ 
changed. In parenthesis, the values related to the cubic 
lattice constant are presented to make comparison easier. 
In the AFM structure, these values are smaller by about 
0.05-0.1% than the cubic lattice constant. In the AFM-II 
structure, the a lattice constant elongates by 0.3% and 
c/2 shortens by 0.9% in comparison to the cubic struc¬ 
ture. 

Since the interatomic distances and the magnetic mo¬ 
ments stay almost unchanged, also the interatomic forces 
and the PDOS presented in Fig. [6] change very weakly 
(apart from small differences observed around 4.7 THz). 
Both of them are similar to the spectrum calculated for 
the FM phase and the shape of the partial Fe PDOS is 
also the same. The small shift to lower frequencies is 
caused by lattice constants used in the calculations of 
AFM-I and AFM-II, which differ slightly from experi¬ 
mental parameters applied for the cubic structure. 

In the present studies, no stresses in the structure are 
assumed, in contradiction with a previous publication 
where authors speculate about the existence of AFM-I 
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FIG. 5. (Color online) The phonon density of states calculated for the LI 2 structure of FesPt and FePt 3 . 


and AFM-II ordering under the uniaxial and tetragonal 
stresses.^ The central issue of that first principle calcu¬ 
lations was to investigate the influence of strain on the 
magnetic state. 

C. Low-temperature phase of FegPt 

The structure of the low-temperature phase of stoi¬ 
chiometric FeaPt ordered alloy has not yet been exactly 
described. X-ray diffraction investigations performed by 
different groups on polycrystals with similar composition 
and degree of atomic order so far did not yield conclusive 
results on the low temperature structure ! 10 This material 
cannot be investigated by means of neutron diffraction 
since Fe and Pt have almost the same coherent scat¬ 
tering lengths, 9.45 frn and 9.60 fm, respectively. Up 
to now, three types of ma rtensi tic phases have been re¬ 
ported: bcc, bet, and fet ! 9 * 31 * 32 -! The bcc martensite is 
formed when the degree of order of the parent phase is 
very low. If the degree of order is in the intermediate 
range, a thermoelastic transformation from the parent 
phase occurs and the bet martensite is formed. For an or¬ 
dering degree of 0.80, the fet martensite forms through a 
second order transformation. The latest X-ray diffraction 
measurements revealed that the highly ordered FeaPt al¬ 
loy exhibits a martensitic transformation below 60 K and 


the tetragonality parameter c/a of the low temperature 
phase is larger (1.005) or lower (0.94) than 1 for the de¬ 
gree of order 0.88 or 0.75, respectively.^ 

In the present calculations, the deformation caused by 
uniaxial stress leading to the tetragonal fet ( PA/mmm ) 
phase with c/a larger and lower than 1 has been con¬ 
sidered. The optimization of those structures allowed to 
find structures with energy slightly lower than that for 
the cubic phase (see Table IV). The optimization proce¬ 
dure finished when the magnitudes of the H-F forces were 
less than 10 -5 eV/A and the external pressure reached 
zero. After optimization, small stresses are observed in 
both tetragonal structures: 3.83 kb in [100] and [010] 
directions and -6.94 kb in [001] direction for tetragonal¬ 
ity c/a >1 and -0.95 kb and 1.91 kb for c/a <1. The 
tetragonal deformation slightly modifies the magnitude 
of magnetic moments of Fe atoms, mainly those lying in 
the (110) plane. They increase with increasing lattice 
constant c. This kind of tetragonal deformation cannot 
remove the imaginary frequencies (2.552i THz for cubic 
structure) from dispersion curves at the M point, it lays 
even deeper at 2.806i THz. 

To stabilize the crystal structure, we have distorted the 
cubic lattice using the soft-mode polarization vector and 
then optimized the lattice parameters with new symme¬ 
try constraints. Previously, the calculations with the dis¬ 
torted FeaPt crystal have been performed to investigate 
















TABLE III. The space group, nonequivalent atomic positions, magnitude of magnetic moments, and energy per 4 atoms unit 
cell calculated for the ordered FePt 3 alloy with different magnetic orders. 


magnetic 

space 

lattice 

non-equivalent 

magnetic 

energy per 

order 

group 

parameters (A) 

atoms 

moments (ps) 

unit cell (eV) 

FM 

Pm3m 

a = 3.914 

Fe(0,0,0) 

3.258 

0 




Pt(f,f,0) 

0.359 


AFM 

Pm3m 

a = 3.891 

Fe(0,0,0) 

3.325 

0 




Pt(i,|,0) 

0.006 


NM 

Pm3m 

a = 3.887 

Fe(0,0,0) 

0.0 

1.101 




Pt(i,|,0) 

0.0 


AFM-I 

PA/mmm 

a = 5.533 

Fe(0,0,0) 

3.272 

-0.094 



(3.912) 

Fe(|, 1,0) 

3.272 




c = 3.900 

Pt(Ud) 

0.000 


AFM-II 

PA/mmm 

a = 3.926 

Fe(0,0,0) 

3.306 

-0.081 



c = 7.756 

Fe(0,0,|) 

3.306 




(3.878) 

Pt(0,|,|) 

0.000 





Pt(i,§,0) 

0.000 





pta.u) 

-0.080 



Energy (meV) 


0 8 16 24 



FIG. 6. (Color online) The phonon density of state calculated 
for the FePt 3 low-temperature phases. 


the changes in electronic structure induced by the soft- 


mod^ 19 . The resulting crystal structure has a lower total 
energy and reduced symmetry to the tetragonal PA/mbm 
space group (No 127) with a doubled unit cell. The dis¬ 
persion curves calculated using this space group are pre¬ 
sented in Fig. [7] The high symmetry directions of the 
cubic phase (Fig{2j) are used since the geometries of the 
cubic and tetragonal phases are comparable. The number 
of curves is 24 as there are 8 atoms in the primitive cell. 
The experimental data measured at room temperature 
for cubic FeaPt are compatible with our calculated dis¬ 
persion curves. Especially, the TAi mode which is soft 
in the cubic phase is perfectly reproduced by the low¬ 
est acoustic branch calculated for the tetragonal phase. 
Additionally, the measured optic modes with too high 
frequencies in the cubic structure can be pinned to some 
curves obtained from these calculations. 

r M t (x c ) r t (M c ) x t r r 


FIG. 7. (Color online) The phonon dispersion curves calcu¬ 
lated for tetragonal structure of FesPt. The experimental 
data are taken from Ref. 8|. 
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TABLE IV. The comparison of high and low-temperature phases of completely ordered FesPt. 


magnetic 

order 

space 

group 

lattice 

parameters (A) 

non-equivalent 

atoms 

magnetic 
moments (^b) 

energy per 
unit cell (eV) 

above M s 

Pm3m 

a = 3.737 

Fe(0,0,0) 

Pt(i,i,0) 

2.710-2.714 

0.349 

0 

below M s 

P4/mbm 

a = 5.253 

Pt(0,0,0) 

0.346 

-0.019 



(3.910) 

Fe(0,|,0) 

2.704 




c = 3.789 

Fe(0.2275,0.2725,0) 

2.721-2.728 



P4/mmm 

a = 3.701 

Fe(i,0,|) 

2.711 

-0.008 



c = 3.813 

Fe(0,§,§) 

2.710 





Fe(|,|,0) 

2.734 





Pt(0,0,0) 

0.357 



P4/mmm 

a = 3.776 

Fe(|,0,0) 

2.729 

-0.016 



c = 3.664 

Fe(0,§,±) 

2.729 





F e (|, |,0) 

2.669 





Pt(0,0,0) 

0.352 



The anomalous temperature dependence of the PDOS 
in the FeaPt ordered alloys has been studied using the 
NIS technique.^ The presence of the soft mode leads 
to enhancement of the intensities of the low-frequency 
vibrations so in the low-frequency region the PDOS pro¬ 
ceeds to zero not in a parabolic but linear manner. In Fig. 
[8j the Fe PDOS measured by Wiele et al} 1 at T = 295 
and 75 K are presented together with the spectra calcu¬ 
lated for the cubic Pm3m and tetragonal PA/mbm struc¬ 
tures, respectively. The agreement between the theory 
and experiment is much better for the tetragonal phase 
measured at 75 K than for the cubic one at room tem¬ 
perature. The positions and intensities of three peaks 
observed in the experiment correspond very well to those 
found theoretically. Also, the shoulder at about 2 THz 
originating from the soft mode is very well reproduced in 
calculations. It clearly demonstrates that the tetragonal 
structure derived from the soft-mode is an important con¬ 
stituent of the low-temperature phase. It may be related 
to static precursors of the martensitic phase observed by 
X-ray, electron, and neutron diffraction and manifested 
as tetragonally strained structures growing in the austen¬ 
ite phase approaching the structural transformation 3335 !. 


It should be noted that this type of crystal deforma¬ 
tion may be driven by the electron-phonon interaction. 
According to previous theoretical studies on FeaPt, the 
condensation of the soft-mode opens a pseudogap in the 
density of states and leads to the extended reconstruc¬ 
tion of the Fermi surface due to nesting effect^. Such 
electron-phonon coupling involving the soft-mode distor¬ 
tion and changes in valence electron density may be a 
characteristic feature of other Invar alloys, which exhibit 
martensitic transformation at low temperatures. 


Energy (meV) 

0 8 16 24 32 40 48 



FIG. 8. (Color online) The phonon density of states calcu¬ 
lated for the tetragonal P4/mbm structure of FesPt. The 
experimental points are taken from Ref. HD¬ 


TV. CONCLUSIONS 

To complete and systematize the understanding of the 
ordered FeaPt, FePt, and FePt 3 alloys, first-principles 
studies of their structural, magnetic, and dynamical 
properties have been performed. The obtained results 
have been compared with existing experimental data sup¬ 
plemented by the new partial Fe PDOS measured in thin 
Fe-Pt films of Ll 2 structure using the nuclear inelastic 
scattering technique. The calculations involve the room 
temperature as well as the low temperature phases of 
Fe-Pt alloys. 

Conclusions about the room temperature structures of 
ordered Fe-Pt alloys can be summed up as follows: 
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(i) The phonon dispersion relations calculated for FePt 
and FePt 3 show very good agreement with inelastic neu¬ 
tron scattering data measured at room temperature. The 
paramagnetic FePt 3 systems were modelled using static 
FM and AFM arrangements of local magnetic moments 
on Fe atoms and both approaches reproduce the exper¬ 
imental data much better than the calculations for the 
non-magnetic state. In the cubic Ll 2 structure of FesPt, 
the observed TAi soft mode at the M point of the Bril- 
louin zone leads to some discrepancy between the calcu¬ 
lated and measured phonon branches. 

(ii) On the basis of the analysis of different variants of 
calculated partial PDOSs, the experimentally unobserv¬ 
able anisotropy of some sites or atoms is noticed. The 
anisotropy of the partial Fe and Pt PDOSs appears ev¬ 
idently in the tetragonal FePt, however, it is also ob¬ 
served for the atomic sites of the tetragonal symmetry in 
the cubic FesPt and FePt 3 . For these cubic structures 
the PDOS asymmetry, which is not important for a bulk 
material, can influence the properties of nanocrystals or 
surfaces. 

(iii) The Fe projected PDOS calculated for the FesPt 
and FePt 3 cubic structures were compared with new nu¬ 
clear inelastic scattering measurements for thin films at 
room temperature. For FePt 3 , we found very good agree¬ 
ment between experiment and theory although additional 
broadening of the measured PDOS is observed. The 
soft-mode behavior observed in FesPt generates some 
disagreement between the theoretical and experimental 
PDOS, mainly in the low-frequency region. 

The chemical ordering in conventional bulk Fe-Pt al¬ 
loys becomes very slow at low temperatures making 
the experimental studies of the low temperature struc¬ 
tures very difficult. Nevertheless, this region was ex¬ 
tensively investigated for the Fe-rich compositions, as 
these alloys show the Invar effect and undergo a marten¬ 
sitic transformation, and for the Pt-rich compositions, 
where the paramagnetic-antiferromagnetic phase transi¬ 
tion is observed. The results of our research on the low- 
temperature phases of FesPt and FePt 3 can be summa¬ 
rized as follows: 

(i) In FePt 3 , the crystal structure parameters, mag¬ 
netic moments, and total energies were determined for 
two antiferromagnetic arrangements that induce the 
tetragonal distortion (P4/mmm). The ground state en¬ 
ergies of AFM-I and AFM-II configurations are nearly 


the same and slightly lower than the energies of the hy¬ 
pothetical FM and AFM phases. The magnitude of mag¬ 
netic moments calculated on Fe atoms is about 3.3 /j,b for 
AFM-I and AFM-II. In contrast, the measured magnetic 
moment of AFM-II phase is 2 /r# and is equal to the value 
in disordered Fe-Pt alloys. It proves the previous findings 
describing the AFM-II phase observed in FePt 3 as formed 
by antiphase domainsPOur calculations showed that any 
antiferromagnetic structure is energetically favoured and 
the changes induced by these two magnetic configura¬ 
tions affect very weakly the lattice dynamical properties 
leading to very similar phonon spectra. 

(ii) In FesPt, the low-symmetry tetragonal structure 
( P^/mbm ) has been derived using the soft-mode polar¬ 
ization vector. The energy of this structure is lower than 
the energies of two tetragonal fct phases (PA/mmm) ob¬ 
tained by uniaxial stress deformation. A few meV per 
unit cell differences between them seems to be irrelevant 
for the stabilization of the low temperature phase, how¬ 
ever, it was shown that the dispersion curves and the 
Fe projected PDOS calculated for the PA/vibm symme¬ 
try agree very well with the inelastic neutron scattering 
and nuclear inelastic scattering (at T = 75 K) data, re¬ 
spectively. It shows that although the martensitic phase 
transition is not a group-subgroup transformation, the 
tetragonal structure derived from the soft mode plays 
an important role in this transition as the main con¬ 
stituent of the low-temperature phase of FesPt. This 
is also the case in the martensitic transformation of NiTi 
alloys where the soft mode of the austenite phase leads to 
either an intermediate incommensurate phase locked into 
a trigonal R phase or to an orthorhombic phase, which 
in turn, creates a low-frequency mode and favours the 
monoclinic martensitic phase 
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